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ABSTRACT 


Hydrodynamic  loads  induced  on  large  underwater  structures 
by  Impulsive  motion  may  be  significant  design  factors.   Such 
loads  may  be  induced  by  earthquake  excitation  or  may  result 
from  acceleration  produced  while  lifting  an  object  in  the 
sea. 

The  theoretical  approach  to  the  calculation  of  these 
loads  is  outlined  and  numerical  results  are  presented  for 
several  submerged  configurations  including  a  vertical  cone, 
a  sphere,  and  a  vertical  circular  cylinder. 

Numerical  results  for  these  submerged  structures  are 
presented  in  the  form  of  a  dimensionless  load  parameter  or 
added  mass  coefficient.   Results  corresponding  to  a  number 
of  different  water  depths  are  presented  to  show  the  rather 
sizable  effect  of  the  relative  water  depth  on  the  hydro- 
dynamic  force.   It  is  shown  that  the  effect  of  the  free 
water  surface  is  to  reduce  the  hydrodynamic  loads  in 
comparison  to  the  corresponding  infinite  depth  values. 

Experimental  results  obtained  by  vibration  testing  are 
presented  for  a  submerged  sphere,  a  cone  and  a  vertical 
cylinder.   These  results  show  excellent  agreement  with  the 
theoretical  results. 
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NOMENCLATURE 

Symbol  Description 

a  =  characteristic  dimension  of  the  body 

A..  =  added  mass  coefficient 

d  =  depth  of  submergence 

f(5n y Zo> ?o)  =  source  strength  function 

f  =  force  vector 

g  =  acceleration  of  gravity 

G  =  Green's  Function 

h  =  depth  of  water 

K  =  modified  Bessel  function  of  the  second  kind 

n  of  order  n 

T  .=  distance  of  line  of  action  above  bottom 

$1  =  moment  vector 

n  =  unit  normal  vector 

p  =  dynamic  pressure,  dimensionless 

P  =  dynamic  pressure 

->■ 

q  =  fluid  velocity  vector 

S  =  surface  of  body 

t  =  time 

U  =  velocity  vector  of  object 

x,,Xp,x.~  =  coordinates,  dimensionless 

2-, 

u  =a  a/g 

C-i  s£p>5o  =  coordinates  of  a  point  on  the  surface  of 

-*  the  object,  dimensionless 

p  =  fluid  density 


Symbol  Description 

a  =  frequency 

<}>  =   velocity  potential,  dimensionless 

*  =  velocity  potential 

ift  =  angular  vector  velocity  of  body 
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I.   INTRODUCTION 

When  an  object  is  accelerated  through  a  fluid  there 
are,  in  general,  two  types  of  forces  that  are  recognized, 
one  being  a  drag  component  and  a  second,  the  inert ial 
component.   The  nature  of  the  flow  produced  by  the  time 
dependent  motion  of  an  immersed  rigid  object  and  the  relative 
contribution  of  these  two  components  of  force  is  generally 
considered  to  be  strongly  dependent  on  the  amplitude  of  the 
relative  fluid  motion  in  comparison  to  the  characteristic 
lineal  dimension  of  the  object.   For  example,  Keulegan  and 
Carpenter  [Ref.  1]  found  that,  for  harmonic  motion  of  the 
fluid  past  a. fixed  circular  cylinder,  major  flow  separation 
did  not  occur  and  the  forces  were  well  represented  by 
potential  flow  values  provided  the  amplitude  of  the  motion 
was  less  than  about  a  half  diameter.   If  the  amplitude 
is  small  enough  the  fluid  does  not  move  in  one  direction 
far  enough  for  separation  and  wake  development  to  occur 
and,  accordingly,  an  inviscid  flow  analysis  is  a  valid 
mathematical  model. 

Accelerated  motion  of  a  solid  body  in  a  fluid  occurs 
for  example  when  a  bottom  mounted  structure  is  excited  by 
an  earthquake.   Other  applications  include  the  accelera- 
tion of  ships  and  accelerated  motions  of  objects  being 
lifted  in  the  ocean  environment.   Southerland  [Ref.  2]  has 
described  the  effect  of  added  mass  on  the  mechanics  of 
lifting  such  loads  as  a  submarine  rescue  chamber. 
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John  [Ref.  3]  has  shown  that  the  velocity  potential 
associated  with  the  harmonic  motion  of  a  rigid  object 
immersed  in  a  fluid  with  a  free  surface  may  be  represented 
by  a  surface  distribution  of  point  wave  sources  for  arbi- 
trary values  of  frequency.   The  source  distribution  (or 
weighting)  function  is  obtained  in  this  formulation  from 
the  solution  of  an  integral  equation  resulting  from  the 
kinematic  boundary  condition  applied  on  the  immersed  surface. 
Using  this  technique,  numerical  results  have  been  obtained, 
for  example,  for  a  semiellipsoid  oscillating  in  a  free  sur- 
face in  water  of  infinite  depth  by  Kim  [Ref.  4],  the 
intended  application  being  to  surface  ships.   Numerical 
results  for  a  bottom  mounted  hemisphere  have  been  presented 
by  Garrison  and  Seetharama  Rao  [Ref.  5]  and  results  for 
several  axisymmetric  configurations  by  Milgram  and  Halkyard 
[Ref.  6].   Garrison  and  Chow  [Ref.  7]  have  recently  presented 
both  theoretical  and  experimental  results  for  wave  forces 
acting  on  a  fixed  submerged  oil  storage  tank. 

In  the  work  cited,  the  hydrodynamic  force  coefficients 
are  represented  as  a  function  of  a  dimensionless  frequency 
parameter  (or  equivalent ly,  a  dimensionless  wave  length 

parameter)  of  the  form  v  =  o  a    where  a   denotes  the 

_  S 
frequency  of  the  motion,  a  denotes  the  characteristic  lineal 

dimension  of  the  submerged  object  and  g  denotes  the  gravi- 
tational constant.   This  dimensionless  parameter  is  essentially 
a  form  of  the  Froude  number.   In  all  cases,  however,  the 
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theory  has  a  fundamental  limitation.   The  source  potential 

oscillates  with  a  wave  length  inversely  proportional  to 

the  dimensionless  frequency  of  the  motion  v,  and  consequently, 

a  reasonable  partition  size  of  the  immersed  surface  limits 

the  validity  of  the  numerical  scheme  to  small  to  moderate 

values  of  frequency.   Specifically,  it  has  been  found 

using  the  Haskind's  relations  [Ref.  8]  and  an  energy  balance 

for  purposes  of  checking  accuracy,  that  results  can  be 

obtained  with  reasonable  accuracy  up  to  about  v  =  2.0. 

Beyond  this,  the  accuracy  of  the  numerical  results  becomes 

questionable . 

Although  the  range  v  =  0  to  about  2.0  includes  most 
applications  in  ship  hydrodynamics  and  wave  forces  acting 
on  large  structures,  the  frequencies  encountered,  for 
example,  in  earthquake  excitation  corresponds  to  rather 
large  values  of  v  in  the  case  of  large  structures.   Accord- 
ingly, the  development  of  an  alternate  method  which  includes 
both  the  effect  of  the  free  surface  and  bottom  on  the 
hydrodynamic  forces  and  is  valid  for  large  values  of  the 
dimensionless  frequency  v  is  the  primary  objective  of  the 
present  thesis.   Therefore  an  attempt  has  been  made  herein 
to  obtain  numerical  results  for  several  different  objects 
for  the  asymptotic  case  of  very  large  values  of  v,  i.e., 
for  the  infinite  Froude  number  case. 

High  frequency  added  mass  coefficients  for  two-dimensional 
shapes  oscillating  on  a  free  surface  have  been  presented 
in  a  series  of  papers  by  Landweber  et  al .  [Ref.  9,  10,   1, 
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and  12].   The  motivation  for  this  work  was  to  provide 
hydrodynamic  coefficients  for  two-dimensional  ship  forms. 
For  three-dimensional  objects  the  only  results  known 
are  those  of  Waugh  and  Ellis  [Ref.  13]  for  the  case  of  a 
sphere  in  infinite  depth  water  and  Jacobsen  [Ref.  14] 
for  the  case  of  the  cylindrical  pier  in  various  water 
depths.   In  the  former  case,  the  theoretical  method  consisted 
of  using  successive  image  doublets  to  account  for  the  effect 
of ' the  free  surface.   However,  the  case  of  vertical  motion 
only  was  studied  and  their  method  was  limited  to  a  sphere. 
In  the  latter  case,  the  solution  was  represented  by  a  series 
of  Bessel  functions  and  is  applicable  to  a  vertical  circular 
cylinder,  extending  between  the  bottom  and  free  surface. 
However,  the  numerical  computations  were  made  using  tables 
of  Bessel  functions  with  complex  arguments  which  have 
proven  misleading.   In  fact,  Jacobsen' s  numerical  results 
are  in  error  in  several  instances  and  therefore,  his 
analysis  as  well  as  the  corrected  formulae  are  presented 
herein. 

In  this  thesis  the  theoretical  formulation  of  the 
impulsive  hydrodynamic  forces  acting  on  a  rigid  object  of 
arbitrary  shape  is  presented.   The  effect  of  the  bottom 
and  free  surface  is  taken  into  consideration  and  the  force 
coefficient  is  shown  to  depend  on  a  frequency  of  oscillation 
parameter.   The  problem  is  first  formulated  for  arbitrary 
values  of  the  frequency  and  the  asymptotic  form  of  the 
solution  for  both  large  and  small  values  of  the  frequer  y 
is  discussed.   A  computer  method  for  numerical  evaluation 
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of  the  force  coefficient  for  submerged  objects  of  arbitrary 
shape  is  discussed  and  numerical  results  are  presented 
for  a  number  of  specific  geometric  shapes  including  a 
circular  cylinder,  cone,  and  sphere.   High  frequency 
experimental  results  are  presented  for  a  vertical  circular 
cylinder,  sphere,  and  cone.   These  experimental  results 
are  compared  with  the  theoretical  results. 
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II.   THEORETICAL  CONSIDERATIONS 

A.   BOUNDARY  VALUE  PROBLEM 

The  geometry  involved  in  the  problem  under  consideration 
is  shown  in  Figure  1.   A  rigid  object  of  arbitrary  shape 
having  characteristic  lineal  dimension  a  is  submerged  to  a 
depth  d  in  water  of  depth  h  and  caused  to  move  in  a  small 
amplitude  harmonic  motion.   The  instantaneous  motion  of  the 
rigid  object  is  specified  by  the  lineal  velocity  U  and 
angular  velocity  ft .   The  object  may  intersect  the  bottom 
or  free  surface. 

Assuming  an  unseparated,  incompressible  and  irrotational 

flow,  a  velocity  potential  $(x,t)  defined  by 

-> 

q  =  V$  (1) 

-*- 
where  q  denotes  the  fluid  velocity  vector,  may  be  introduced 

and  must  satisfy  the  continuity  equation, 

V2$  =  0  (2) 

within  the  fluid  region.   The  bars  over  the  spatial  variables 

-> 

denote  dimensional  quantities  and  x  =  ix,  +  jx?  +  kx_. 

On  the  rigid  bottom  surface  described  by  x2  =  -h  the 
fluid  velocity  in  the  vertical  direction  must  vanish  and, 
accordingly,  the  corresponding  kinematic  boundary  condition 
on  the  velocity  potential  is 


3$  (x-,  ,-h,x,,t)  = 

;2 


.1,-n,x3,i:;  =  Q  ^) 

3X, 


15 


The  velocity  of  a  point  on  the  immersed  surface  is  given 

by 

v1  =  8  +  fl  x  r  (4) 

-*■ 
where  r  denotes  the  position  vector  as  shown  in  Pig.  1.   The 

appropriate  kinematic  and  dynamic  boundary  condition  on  the 

immersed  surface,  therefore,  takes  the  form 

*!<*»*>  =  f-n"  (5) 

9n 

where  n  =  in-,  +  jn2  +  kn~  denotes  the  unit  normal  vector  on 
the  surface  of  the  object  and  is  directed  outward  into  the 
fluid  region. 

Assuming  the  amplitude  of  the  motion  of  the  object  to 
be  small,  the  kinematic  and  dynamic  boundary  condition  applied 
at  the  free  surface  may  be  linearized  to  give  the  well-known 
free  surface  boundary  condition: 

2iJx1,0,x3,t)  +  1  3_£(x190,x39t)  =  0         (6) 

3X2  g  3t2 

Furthermore,  if  the  object  oscillates  with  frequency  a,  the 
time  dependence  of  the  velocity  potential  may  be  separated 
and  written  as  Re[<}>  (x,  ,Xp,x_)e~la  ]  where  Re  denotes  the  real 
part.   In  which  case  Eq.  (6)  becomes 

f^(xl'°-x3)  -  ^   ♦(x1>0,x3)  =  0   '         (7) 

where  the  unbarred  quantities  represent  the  coordinates   ade 
dimensionless  with  a. 
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Fig.    1   -  Definition   of   Geometry 
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In  general  it  is  apparent  from  Eq.  (7)  that  the  solution 

to  the  boundary  value  problem  established  must  depend  on  the 

2— 
dimensionless  frequency  parameter  a  a/g.   However,  there  are 

two  extreme  cases  where  the  results  become  independent  of  the 

2— 
frequency.   These  correspond  to  the  case  where  a   a/g  ■+  °°  in 

which  case  Eq.  (7)  becomes 

♦(x1,0,x3)  =  0   •  (8) 

2— 
For  the  other  extreme  where  a  a/g  ■*  0  the  free  surface  appears 

as  a  rigid  boundary  as  Eq.  (7)  becomes 


|1l(x1.o.x3)  .  0  (9) 

d  X/^ 


If  consideration  is  restricted  to  either  of  the  boundary 
conditions  Eq.  (8)  or  Eq.  (9)  the  velocity  potential  may  be 
written  as  the  sum 

4>(x,t)  =  W±(t)  ♦±(x),    i  =  1,2,. ..6        (10) 

where  each  term  in  the  sum  represents  the  potential  associated 

A-  i_ 

with  the  i   component  of  body  motion.   (i  =  1,2,3  refers  to 
linear  motion  components  in  the  x-.,Xp,x~  directions, 
respectively,  and  i  =  4,5,6  refers  to  the  angular  motion 
components  about  the  x-,  ,Xp,x^  axes,  respectively.)   The  time 
dependent  function  W  is  defined  as 

W, (t)  =  U, (t)  a 

Wi+3(t)  =  D±(t)  i"2 
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In  which  the  symbols  U.  and  ft.  denote  the  components  of  the 
lineal  and  angular  velocity  vectors  of  the  body,  0  and  Q, 
respectively.   The  potentials  <j>.  are  dimensionless  quantities 

It  is  important  to  note  that  the  application  of  either  of 
the  two  boundary  conditions,  Eqs.  (8)  and  (9)s  makes  the 
problem  homogeneous  in  time.   Consequently,  the  assumption 
of  harmonic  motion  may  be  relaxed  and  the  functions  W. (t )  may 
be  allowed  to  be  arbitrary  functions. 

Substituting  Eq.  (10)  into  Eqs.  (2),  (3),  (5)  and  (8)  or 
(9)  and  introducing  the  dimensionless  variables  defined  as 

x.  =  x./a,   h  =  h/a,   d  =  d/a  (12) 

the   boundary  value  problem  for   $.    (i   =   1,2, ...6)   may  be 
specified   as    follows: 


V   <{.jL(x1,x2,x3)    =   0 


a<f>. 


<J>i(x1,0,x3)    =   0,    (o-m»)      or     _!(x150,x3)    =   0,    (a-0) 


3<j>.  (x,  ,-h,x,) 

jf-     1  3     =   0  (13  a-d) 

- —  =  m.    (on  the  immersed  surface) 

9n  i 


where  m..    =  n.    ,    i   =   1,2,3 


mi,  =    (xn+d)n^,   -   x.-,n0 
q  2  3  3   2  (14   a-d) 

nir  =   x0nn    -  x,n_ 
o  5   1  13 

m,-  =   x,n2   -    (x2   +  d)n. 
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The  boundary  conditions  (13b)  represent  two  different 
possibilities  on  the  free  surface,  the  first  being  the  high 
frequency  condition  and  the  second  the  rigid  boundary  (or  zero 
frequency)  condition.   In  as  much  as  the  method  of  solution 
involving  either  boundary  condition  is  similar,  both  solutions 
are  developed  simultaneously  in  the  following. 

B.   GREEN'S  FUNCTION  SOLUTION 

The  solution  to  the  boundary  value  problem  (13)  may  be 
carried  out  by  use  of  a  Green's  function.   The  Green's 
function  represents  the  potential  for  a  point  source  and 
these  sources  are  distributed  over  the  immersed  surface 
according  to  the  distribution  function  f(c).   According  to 
this  concept  the  potential  <f>  may  be  written  as 

*i(*}  =  1J7  f  f  V^  G(^)dS  (15) 


y 


s 


where  (O  denotes  a  point  on  the  immersed  surface,  G  denotes 

—  —2 

the  Green's  function  and  dS  =  dS/a  denotes  the  surface  area 

element  made  dimensionless  with  the  characteristic  body 
dimension  a. 

In  order  that  Eq.  (15)  represent  a  solution  to  Eqs.  (13) 
it  is  necessary  that  G  satisfy  the  equation 

v2G(x;t)  =  fi(x-t)  (16) 

where  6(x)  denotes  the  Dirac  delta  function,  as  well  as 
boundary  conditions  (13c)  and  (13b). 
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The  Green's  function  satisfying  Eq.  (16)  as  well  as 
Eq.  (13b)  and  (13c)  is  obtained  by  use  of  successive  images 
of  sources  as  depicted  in  Fig.  2  and  is  given  by 

G(hi)  =  |  -  jJ-  +Z;  [(-Dn  (it-  +  ir->  +  (-i>n+1<ir- +  ir->] 

R    Rl    n=1  R^    R^  R3n    R^ 

(17) 

The  corresponding  form  of  G  satisfying  the  rigid  boundary 
condition  Eq.  (13c)  (o  =  0  on  x,  =  0)  is  also  obtained  by 
successive  images  and  is  similar  to  Eq.  (17)  except  that  all 
of  the  sources  have  positive  strengths.   For  this  case,  the 
Green's  function  takes  the  form 


«iit)  -  |  +  ^  +  ^  <j£-  ♦  ^  ♦  rJ-  +  rM  (18) 

n    Jn     n     n 

where  in  both  Eq.  (17)  and  (18) 


a)   =  y  (x1  -  c-l)2  +  (x3  -  ?3)2 


R   =~Vo)2  +  (x2  -  ?2)2 

R1   =  V"2  +  (x2  +  ^2)2 

R2   =  V  u>2  +  (x2  -  2nh  -  c2)2  (19) 

n 

R~   ="\/w2  +  (x2  +  2nh  +  C2)2 
^n 
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Since  G  satisfies  Eq.  (13a-d)  the  potential  <J>.  as  given  by 
Eq.  (15)  "must  also  satisfy  these  conditions.   Thus,  it  remains 
to  select  the  function  f  such  that  Eq.  (13d)  is  satisfied. 
Application  of  this  boundary  condition  yields  the  following 
Fredholm  integral  equation  of  the  second  kind 

-f±(x)  +  ±  J  j  f±(t)   f£(x;°   dS  =  2m1(x)    ,  i  =  1,2, ...6  (20) 

o 

which  is  to  be  solved  for  f(5).   This  integral  equation  is  to 
be  satisfied  at  all  points  x  on  the  immersed  surface. 

Because  of  the  form  of  G,  Eq.  (20)  is  rather  complex  and 
it  is  therefore  necessary  to  carry  out  the  solution  to  the 
integral  equation  numerically.   Proceeding  in  this  direction 
the  immersed  surface  may  be  partitioned  into  N  area  elements 
of  size  AS.  and  the  integral  in  Eq.  (20)  written  as  the  sum 

-f    +  f    a. .  =  2m    ,     n  =  1,2, ...6  (21) 

n±         n.   ij     n1  *         '  5 

where  i,j  =  1,2, ...N  and  n  refers  to  the  six  possible  modes 
of  motion  of  the  object.   The  transformation  from  the  integral 
equation  to  a  summation  is  possible  since  f  is  a  well-behaved 
function.   It  follows  from  Eq.  (20)  that  the  coefficient  matrix 
is  defined  by 
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=  iff    pi   ;t)  (   } 

•J 

The  function  3G/an  in  Eq.  (22)  is  obtained  in  a  straight- 
forward manner  by  differentiation  of  Eq.  (17)  or  (18). 

Once  the  coefficient  matrix  a,,  is  obtained  from  Eq.  (22) 
a  standard  digital  computer  subroutine  may  be  used  to  invert 
the  matrix  equation,  Eq.  (21),  to  obtain  the  source  strength 
f  at  points  on  the  immersed  surface.   For  purposes  of  evalua- 
tion of  G  and  3G/an  by  use  of  (17)  and  (18)  it  was  found  that 
about  15  terms  were  adequate.   The  data  presented  herein, 
however,  were  obtained  using  30  terms. 

C.   HYDRODYNAMIC  FORCES  AND  MOMENTS 

The  dynamic  fluid  pressure  is  obtained  from  the  linearized 
form  of  Bernoulli's  equation  as 


P--.P&  (23) 


where  the  velocity  squared  term  has  been  neglected. 
Corresponding  to  the  i    component  of  body  motion,  the  pressure 
is  obtained  by  use  of  Eq.  (10)  and  (23)  as 

?±   =   -pw.<}>.  (no  sum)       i  =  1,2, ...6      (24) 

The  pressure  may  further  be  written  in  dimensionless  coeffi- 
cient form  as 

P±(x)  =  — =  -4>.(x)  ,    i  =  1,2, ...6      (25) 

pw. 
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The  hydrodynamic  force  acting  on  the  object  is  obtained 
by  carrying  out  the  integration  of  the  pressure  over  the 
immersed  surface.   Using  Eq.  (24)  the  force  vector  associated 

i.  ■u. 

with  the  j    component  of  body  motion  is 

*J  =  -^2  j  J   Pj-ndS,      j  =  1,2,... 6         (26) 


The  corresponding  result  for  the  moment  vector  is 

P.  r  x  n  dS,   j  =  1,2,. ..6      (27) 


-3  " 

M.  =  -a-3 


Using  Eq.  (24)  as  well  as  the  definition  of  m.  as  given  in 
Eq.  (14),  the  i   component  of  force  associated  with  the 
j    component  of  acceleration  is  given  in  dimensionless 
coefficient  form  as 


Fii_  rn 


13    -SI 


4>,m.  dS,      i  "  1»2»3      (28) 
31  j  =  1,2, ...6 


and  the  corresponding  expression  for  moment  is 


where  w  is  defined  by  Eq .  (11). 

Equations  (28)  and  (29)  define  the  inertia  (or  added  mass) 
tensor  having,  in  general,  36  elements.   The  dimensionless 
coefficients  A. .  defined  by  Eq.  (28)  represent  the  dimension- 
less force  coefficient  and  Eq.  (29)  represents  the  dimension- 
less moment  coefficient.   These  coefficients  are  generally 
referred  to  as  added  mass  (or  moment  of  inertia)  coeff   ents. 
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Although  these  coefficients  are  made  dimensionless  by  use  of 
a  in  the  definitions  (28)  and  (29),  there  are  cases,  depending 
on  the  body  geometry,  where  other  length  scales  are  more 
appropriate.   It  is  also  standard  practice  in  many  cases  to 
use  the  displaced  fluid  volume.   However,  in  all  numerical 
results  presented  herein  the  definition  of  A. .  is  specified 
on  the  figure  so  that  although  the  definitions  vary,  no 
confusion  should  result. 

For  purposes  of  numerical  evaluation  of  the  coefficients 
given  in  Eqs .  (28)  and  (29)  the  integrals  may  be  written  as 
the  summation,  where  the  integral  sign  has  been  replaced  by 
the  sum  on  the  k  index. 

A..  =  U  m  )   AS    ,      1»J    1,2, ...6       (30) 
10      J  l  *   *        k=l,2,...N 

However,  in  order  to  apply  Eq.  (30)  it  is  necessary  to  first 

evaluate  <J> .  at  all  of  the  nodal  points  on  the  immersed 
j 

surface.   For  this  purpose  the  surface  integration  in  Eq.  (15) 
is  also  converted  to  the  summation 

(*1}k  =  ek*(fA   '      j  =  1»2"'"6        (3D 

J  K     Kl      J  *  k,£  =  1,2,.  ..N 

where   the   coefficient   matrix   is   defined   as 


6k*   "  TfiT, 


G(x.;?)     dS  (32) 

'AS„ 


Once  the  coefficient  matrix  is  evaluated  by  Eq.  (32), 
the  summation  indicated  in  Eq.  (31)  may  be  carried  out  t 
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determine  <j>  .  at  all  nodal  points  on  the  immersed  surface. 
The  coefficients  A.,  are  then  obtained  by  use  of  Eq.  (30). 
Although  there  are  36  possible  added  mass  coefficients 
defined  by  Eqs.  (28)  and  (29),  many  are  zero  on  account  of 
symmetry  and  most  off-diagonal  elements  are  small.   Also, 
even  for  a  body  of  arbitrary  shape  A.,  is  symmetrical  so  that 

at  most  only  21  of  the  coefficients  need  be  evaluated.   That 

A 
is,  Eq.  (13e)  along  with  Eq.  (28)  and  (29)  gives 


Au 


s  *Jmi  dS   '  I  |c   *J  5TT  dS        <33) 


By  use  of  the  Green's  theorem  applied  to  the  fluid  region 
external  to  the  body,  it  can  be  shown  that 

since  on  the  free  surface  either  <J>  =  0  or  3<(>/3n  =  0  (depending 
on  the  boundary  condition  being  enforced)  and  3<(>/3n  =0  on 
the  bottom.   Therefore,  in  view  of  Eq.  (3*0  and  (33)  it  is 
apparent  that  the  added  mass  tensor  is  symmetrical, 


Au  -  V  (35) 


D.   JACOBSEN'S  SOLUTION 

At  this  juncture  it  is  appropriate  to  comment  on  the 
analytical  solution  presented  by  Jacobsen  [Ref.  2]  for  the 
special  case  of  a  vertical  circular 'cylinder  extending  between 
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a  rigid  bottom  and  a  free  surface.   This  solution  represents 
the  only  known  analytical  result  for  the  problem  under 
consideration  and,  therefore,  is  of  considerable  value  as  a 
point  of  comparison  for  the  numerical  results  presented 
herein.   However,  it  was  found  that  a  number  of  errors  were 
present  in  Jacobsen's  work  and  it  was  necessary  to  re-develop 
and  re-evaluate  the  equations.   The  expression  for  the  added 
mass  coefficient  associated  with  horizontal  acceleration  in 
terms  of  the  notation  of  the  present  paper  and  in  terms  of 
Bessel  functions  having  real  as  opposed  to  imaginary  argu- 
ments appears  as 


-B  K.(S£0  (36) 

n   n   1  2h 


where  A,-,  denotes  the  added  mass  coefficient  made  dimension- 
less  with  the  displaced  volume,  the  factor  B  is  defined  as 

g     =  l/n (37) 

n        „    /nirx        2h  ,nn » 

Ko(2h}    +  nT  Kl    (2h} 


and  K  and  K,  denote  modified  Bessel  functions  of  the  second 
o      1 

kind  of  order  zero  and  one,  respectively.   This  expression 
agrees  with  Jacobsen's  corresponding  expression  with  the 
exception  of  a  factor  of  2.0  and  the  type  of  Bessel  functions 
involved.   It  was  found  that  more  serious  sign  errors  occurred 
in  Jacobsen's  expression  for  the  line  of  action  of  the  hori- 
zontal hydrodynamic  force  and,  accordingly,  his  numerical 


All = 

l6h 
M3 

>■• 

j 

TT 

n=l 

,3,5 
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results  are  also  in  error.   The  correct  expression  for  the 
elevation  of  the  line  of  action  (or  center  of  added  mass)  in 
terms  of  the  notation  of  the  present  paper  is 


i       n+1  o 

J 


i   =  n=1>3'l> (38) 

/     n  n  lv2h' 


n=l,3,5 

where  i   =  T/h3  T  being  the  elevation  of  the  line  of  action 
above  the  bottom  and  h  the  water  depth. 

For  completeness  the  expression  for  the  dimensionless 
pressure  on  the  surface  of  the  cylinder  is  given  in  terms  of 
the  notation  of  the  present  paper  as 

p  =  _  8h  Cos  9   Z_j   (-1)  2  cos[|  (h+y)]  B  K  (§£)     (39) 
tt        n=l,3,5 

where  6  denotes  the  angle  measured  with  respect  to  the  x-axis 
in  the  x-z  plane. 

Equations  (36)  and  (37)  were  programmed  for  numerical 
evaluation  on  a  digital  computer  so  that  an  adequate  number 
of  terms  could  be  evaluated  without  concern  over  numerical 
error.   These  results  are  plotted  in  Figs.  8  and  9  for 
comparison  with  the  source  distribution  results  developed 
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herein.   It  was  noted  that  numerical  results  obtained  from 
figures  presented  by  Jacobsen  [Ref.  2]  are  in  error  when 
compared  to  numerical  results  obtained  from  Eqs.  (36)  and 
(38). 
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III.   EXPERIMENTAL  APPROACH 

Experiments  were  conducted  to  determine  the  added  mass 
coefficient  of  a  circular  cylinder,  a  cone  and  a  sphere  in 
proximity  of  a  free  surface  and  rigid  bottom.   The  basic 
experimental  approach  consisted  of  mounting  the  test 
object  on  flexible  beams  and  measuring  its  natural  frequency 
of  vibration  in  air  and  water.   The  reduction  in  natural 
frequency  when  placed  in  water  was  then  related  to  the 
added  mass  of  the  water. 

The  test  rig  was  used  to  determine  the  added  mass 
coefficient  corresponding  to  the  horizontal  motion  is  shown 
in  Fig.  3  and  4.   The  sphere  was  mounted  in  this  test  rig 
on  a  thin  6  inch  long  by  2  inch  chord  vertical  strut.   Two 
large  steel  channels  were  placed  perpendicular  to  each 
other  and  welded  to  a  1/2-inch  steel  plate  to  form  a  rigid 
support  for  the  spring-mass  system.   These  channels  were 
supported  at  four  points  on  the  top  edge  of  a  shallow,  water- 
tight box  (24"  x  70"  x  48").   The  cylinder,  cone,  and  sphere, 
mounted  on  a  pair  of  flexible  beams,  were  bolted  to  this 
supporting  structure  as  shown  in  Fig.  3.  (The  cylinder 
is  shown. ) 

Four  strain  gages  were  mounted  on  the  upper  end  of 
one  of  the  flexible  beams  and  connected  into  a  Wheatstone 
bridge.   The  bridge  was  then  connected  to  a  carrier  pre- 
amplifier and  the  frequency  of  the  bridge  output  was  measured 
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Fig.  3  -  Experimental  equipment. 
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by  use  of  a  Monsanto  programmable  counter-timer  and  in  some 
cases  by  the  use  of  a  Visicorder.   A  ten  second  duration 
was  set  on  the  counter-timer  and,  depending  on  the  depth  of 
water  and  test  object,  the  readings  ranged  from  about  110 
to  180  cycles.   The  time  duration  programmed  on  the  counter 
was  highly  accurate  so  that  the  major  error  in  the  process 
was  considered  to  result  from  reading  the  number  of  cycles 
occurring  in  the  10  second  interval  to  integer  values. 
However,  since  the  starting  of  the  interval  was  random, 
this  error  was  minimized  by  repeating  each  test  run  six 
times  and  recording  the  average  reading.   It  was  noted  that 
the  spread  in  readings  from  the  counter  was  never  greater 
than  one  cycle  in  the  ten  second  interval  for  a  given 
configuration. 

A  calibration  curve  for  the  system  was  developed  by 
plucking  the  mass  (test  object)  in  the  air  with  the 
addition  of  pre-measured  weights  attached  to  the  test 
object.   As  a  result  it  was  unnecessary  to  rely  on  the 
theoretical  natural  frequency  equation  for  the  spring- 
mass  system.   In  general,  the  test  procedure  was  considered 
to  be  highly  accurate  and  good  repeatability  was  obtained 
in  all  cases.   The  calibration  curves  are  shown  in  Appendix 
C. 

The  same  procedure  was  used  for  the  determination  of  the 
added  mass  coefficients  for  the  sphere  in  vertical  motion 
but  a  different  apparatus  and  tank  was  employed.   The 
parallel  beam  mount  which  carried  the  rod  with  attached 
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inch  diameter  plexiglass  sphere  is  shown  in  Fig.  5  and 

6.   The  tank  used  in  this  experiment  was  4.0  ft.  deep  and 

represented,  for  practical  purposes,  an  infinite  depth 

case.   The  sphere  was  oscillated  in  the  vertical  direction. 
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IV.   DISCUSSION  OF  RESULTS 

A.   CIRCULAR  CYLINDER 

Theoretical  and  experimental  values  of  added  mass 
coefficients  have  been  generated  for  a  number  of  configurations 
and  are  plotted  in  Figs.  7-20  and  tabulated  in  Tables  1- 
18.   These  include  results  for  a  vertical  circular  cylinder, 
a  sphere,  a  hemispherical  bottom-mounted  vessel,  and  a 
conically-shaped  object. 

The  theoretical  and  experimental  results  for  the  added 
mass  coefficient  for  a  vertical  circular  cylinder  are  shown 
in  Figs.  7-12  and  Tables  1-7.   The  calibration  curves  are 
presented  in  Fig.  7.   Two  different  cylinders  of  the  same 
overall  dimensions  and  materials  were  used.   They  produced 
calibration  curves  which  were  slightly  different -results 
since  their  natural  frequencies  of  vibration  were  different 
and  their  positions  in  the  test  rig  could  not  be  duplicated 
exactly. 

The  added  mass  coefficient  corresponding  to  the  case 
of  the  horizontal  acceleration  as  based  on  the  cylinder 
radius  cubed  is  presented  in  Fig.  8,  both  with  a  free  sur- 
face and  with  a  rigid  upper  boundary.   Results  for  two 
different  grid  sizes  are  shown,  120  and  192  effective  nodal 
points,  and  little  difference  was  found  to  exist.   It  can 
be  seen  that  the  computer  results  for  the  "rigid  boundary" 
case  closely  follow  the  straight  line  relationship  (nh) 
representing  the  classical  two-dimensional  result.   The 
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results  corresponding  to  the  "free  surface"  case  show  that 
the  effect  of  the  free  surface  is  to  reduce  the  added  mass 
coefficient  and  for  larger  values  of  h  the  results  appear 
to  be  closely  approximated  by  the  straight  line  relation- 
ship A1]L  =  Uh  -  3.0). 

The  same  type  of  result  in  a  slightly  different  form  is 
presented  in  Pig.  9.  In  this  figure  the  added  mass  coeffi- 
cient made  dimensionless  in  the  more  common  manner,  by  use 

— 2— 
of  the  displace  volume  ira  h  ,  is  plotted  for  two  different 

grid  sizes  along  with  Jacobsen's  results,  Eq.  (36).   Also, 

the  experimental  results  obtained  from  vibration  testing 

are  presented  for  comparison. 

The  results  show  that  the  effect  of  the  free  surface 
is  to  reduce  the  added  mass  coefficients  and  agreement 
between  the  experimental  and  theoretical  results  is.-  shown 
to  be  excellent.   On  the  same  figure  the  computer  results 
for  the  rigid  boundary  case  are  presented.   The  classical 
result  for  two-dimensional  flow  gives  A, ..  =  1.0  and, 
therefore,  the  deviation  from  this  value  gives  the  percentage 
numerical  error  in  the  computer  results  directly.   The 
figure  indicates  a  maximum  of  about  3.0  per  cent  error  at 
the  extreme  low  values  of  h. 

There  is  no  doubt  that  accuracy  of  the  computer  results 
could  be  improved  by  means  other  than  the  most  obvious: 
increasing  the  number  of  area  elements.   However,  any  gain 
in  precision  is  generally  offset  by  increased  computer 
storage  requirements,  complexity  and  computer  run  tir     The 
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present  method  is  considered  to  be  quite  simple  and  versa- 
tile with  respect  to  configuration  and  yet  possesses  adequate 
accuracy  for  most  purposes. 

For  the  case  of  horizontal  motion,  the  line  of  action 
of  the  horizontal  force  is  also  a  parameter  of  interest. 
When  the  upper  boundary  acts  as  a  "wall",  the  force  acts 
at  the  center,  i.e.,  at  %   *  .5.   However,  for  the  free 
surface  case  the  pressures  are  reduced  near  the  free  surface 
and  the  line  of  action  is  somewhat  below  center.   This  result 
is  presented  in  Fig.  10  as  obtained  using  252  nodal  points 
and  is  compared  with  Jacobsen's  corrected  result,  Eq.  (38). 
The  figure  shows  the  agreement  to  be  excellent. 

A  second  degree  of  freedom  is  also  possible  for  the 
vertical  cylinder.   The  cylinder  may  be  allowed  to  rotate 
about  some  convenient  axis  with  angular  acceleration.   In 
which  case  an  added  moment  of  inertia  may  be  defined. 
Numerical  results  have  been  generated,  using  192  nodal  points, 
for  this  coefficient  for  rotation  about  an  axis  passing 
through  the  center  of  the  cylinder  at  its  base.   These 
theoretical  results  are  presented  in  Fig.  11.   The  results 
show  that  the  free  surface  has  the  effect  of  reducing  the 
hydrodynamic  reaction  when  compared  to  the  rigid  boundary 
case.   This  is  the  expected  result  since  the  free  surface 
tends  to  reduce  the  pressure  differential  near  the  top  of 
the  cylinder. 

Finally ,  in  order  to  show  how  the  maximum  pressure  on 
the  surface  of  the  cylinder  varies  with  depth,  Fig.  1°   as 
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Fig. • 11   -  Added   moment    of   inertia   coefficient 
for   a   vertical   circular   cylinder. 
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vertical  cylinder. 
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been  prepared.   In  this  figure  the  dimensionless  pressure 
coefficient  as  defined  by  Eq.  (25)  is  plotted  for  the  case 
of  horizontal  acceleration.   It  may  be  noted  that  the 
classical  solution  corresponding  to  the  rigid  boundary  case 
yields  the  maximum  pressure  of  p,  =  1.0  over  the  complete 
length.   Fig.  11  shows  that  the  pressure  corresponding  to 
the  free  surface  case  begins  at  zero  and  appears  to  approach 
unity  as  the  distance  from  the  free  surface  and  depth  increase. 
This  is  the  expected  trend  since  the  effect  of  the  free 
surface  diminishes  at  large  distances.   Eq.  (39)  is  plotted 
on  the  same  figure  for  comparison  and  the  agreement  appears 
to  be  quite  good. 

B.   SPHERE 

Results  for  a  sphere  accelerating  both  horizontally  and 
vertically  are  shown  in  Figs.  13-17  and  tables  8-14.   The 
calibration  data  is  presented  in  Fig.  13  (vertical)  and 
Fig.  14  (horizontal).   For  the  case  of  vertical  acceleration 
in  deep  water,  Fig.  15  shows  experimental  results  obtained 
by  vibration  testing  for  comparison  with  theory  and  the 
agreement  appears  to  be  quite  good.   It  may  be  noted  that 
the  effect  of  the  free  surface  is  always  to  decrease  the 
added  mass  coefficients  while  the  proximity  of  a  rigid 
boundary  has  the  opposite  effect.   All  of  the  computer  results 
for  the  sphere  have  been  generated  using  264  nodal  points. 

The  well-known  classical  solution  for  a  sphere  gives  an 
added  mass  coefficient  of  .5  for  the  case  of  a  fluid  of 
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infinite  extent.   Accordingly,  as  d  becomes  large,  the 
coefficient  plotted  in  Fig.  12  is  expected  to  approach 
this  closed  form  result  but  approximately  a  4.0  per  cent 
deviation  is  noted.   This  error  is  due  to  numerical  inaccur- 
acies and  approximations  used  in  the  computer  program. 

The  same  kind  of  results  are  plotted  in  Fig.  16  as 
generated  by  the  computer,  the  only  difference  being  the 
finite  depth,  h  =  4.0.   The  same  general  trends  as  in  Fig. 
12  are  observed  when  the  sphere  is  near  the  free  surface; 
the  rigid  boundary  results  are  greater  than  .5  while  the 
free  surface  results  are  less  than  .5.   However,  as  the  depth 
is  increased  and  the  sphere  approaches  the  bottom,  the 
added  mass  coefficient  increases  in  either  case. 

Theoretical  and  experimental  results  for  the  added 
mass  coefficient  .in  horizontal  acceleration  corresponding 
to  a  bottom  mounted  sphere  is  shown  in  Fig.  17.   The 
experimental  results  were  obtained  by  use  of  the  test  rig 
shown  in  Figs.  3  and  4.   The  sphere  was  mounted  on  the  lower 
plate  by  use  of  a  6.0  inch  strut  and  the  sphere  was  adjusted 
to  within  .02  inches  of  the  tank  floor.   It  may  be  noted 
that  the  results  show  good  agreement  with  theoretical 
results  and  have  the  same  general  trends  as  previous  results 
except  for  the  break  off  in  A,,  for  the  rigid  boundary  case 
as  h  becomes  small.   This  reversal  in  trend  is  caused  by 
the  unwetting  of  the  sphere  occurring  at  values  of  h  less 
than  2.0. 
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Fig.  16  -  Added  mass  coefficient  for  a  sphere 
for  h  =  4.0. 
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C.  HEMISPHERE 

Figure  18  and  Table  15  present  the  added  mass  coeffi- 
cient associated  with  a  bottom  mounted  hemispherical  object 
as  a  function  of  water  depth  for  both  the  rigid  boundary  and 
free  surface  case.   The  results  appear  somewhat  similar 
to  previous  ones  with  the  rigid  boundary  having  the  effect 
of  increasing  the  added  mass  coefficient  and  the  free  surface 
having  the  opposite  effect. 

D.  CONE 

Conically  shaped  bottom  mounted,  surface  piercing 
structures  have  been  proposed  for  use  in  the  Arctic  on 
account  of  the  ice-breaking  capability  of  the  shape.   This 
structure  shown  in  Fig.  19  is  a  conically  shaped  tower 
similar  to  those  which  have  been  proposed  as  drilling  rig 
configurations  for  deployment  in  the  Arctic.   The  computer 
results  were  generated  using  192  effective  nodal  points 
distributed  over  the  immersed  surface.   The  tabular  data 
is  listed  is  Tables  16-18.   The  calibration  curve  is  pre- 
sented as  Fig.  20. 

The  theoretical  and  experimental  added  mass  coefficient 
(based  on  actual  displaced  volume)  for  the  conically  shaped 
tower  is  presented  in  Fig.  19  as  a  function  of  water  depth 
and  these  results  are  in  good  agreement.   The  results  for 
the  rigid  boundary  case  shows  coefficients  which  are  much 
greater  than  those  corresponding  to  the  free  surface  case. 
Moreover,  for  the  free  surface  case  the  line  of  action   f 
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the  hydrodynamic  force  was  found  to  actually  lie  below  the 
bottom.   This  rather  unexpected  result  is  caused  by  the 
fact  that  the  pressures  on  the  lower  portion  contribute  to  a 
moment  opposite  that  induced  by  pressures  on  the  upper  part. 
In  general,  as- the  object  is  accelerated  horizontally, 
positive  pressures  are  produced  on  the  front  side  while 
negative  pressures  occur  on  the  back  side.   Therefore,  since 
the  pressure  as  well  as  area  which  contributes  to  a  negative 
moment  is  much  greater  than  those  tending  to  produce  a 
positive  moment,  the  result  is  a  moment  tending  to  tip 
the  top  of  the  cone  in  the  direction  of  the  acceleration. 
This  moment  is  equivalent  to  applying  the  horizontal  force 
at  a  distance  l   below  the  bottom. 
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V.   CONCLUSIONS 

A  practical  numerical  analysis  has  been  developed  using 
three  dimensional  sources  and  image  sources  to  represent 
the  flow  produced  during  impulsive  motion  of  rigid  struc- 
tures of  arbitrary  shape  immersed  in  water  in  proximity  to 
a  free  surface.   The  bottom  and  free  surface  have  been  taken 
into  account  and  it  is  shown  that  the  proximity  of  the  free 
surface  has  a  sizable  effect  and  always  tends  to  reduce  the 
added  mass  coefficients  in  comparison  to  the  rigid  boundary 
case.   The  experimental  results  obtained  by  vibration 
testing  for  a  sphere,  cylinder,  and  a  cone  were  found  to 
agree  well  with  the  theoretical  results  and  the  results  for 
the  cylinder  are  in  agreement  with  analytical  results. 
The  results  are  directly  applicable  to  such  practical 
problems  as  the  determination  of  earthquake  loading  of  rigid 
submerged  and  semi-submerged  bottom  mounted  structures  as 
well  as  the  calculation  of  required  forces  to  lift  objects 
of  arbitrary  shape. 
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APPENDIX  A 
TABLES  OF  COMPUTER  GENERATED  DATA 

TABLE  1 

CYLINDER:   COMPUTER  GENERATED  DATA  (Horizontal  Motion) 
192  Effective  Points 


h 

A11=P11/pU17Ta2h 

Free  Surface 
Boundary  Cond. 

A11=F11/pU1Tra2h 

Solid 
Boundary  Cond. 

h  An 

Free  Surface 
Boundary  Cond. 

10.0 

0.910 

1.000 

0.472 

9.0 

0.896 

0.999 

0.469 

8.0 

0.879 

0.999 

0.465 

7.0 

0.861 

1.001 

0.461 

6.0 

0.837  • 

1.004 

0.455 

5.0 

0.807 

1.001 

0.449 

4.0 

0.764 

1.014 

0.422 

3.0 

0.699 

1.020 

0.434 

2.0 

0.593 

1.027 

0.425 

1.5 

0.512 

1.031 

0.419 

1.0 

0.402 

1.027 

0.414 

0.5 

0.254 

1.062 

0.412 
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TABLE  2 

CYLINDER:   COMPUTER  GENERATED  DATA  (Horizontal  Motion) 
120  Effective  Points 


h 

— 2~ 
A11=F11/pU1Tra  h 

Free  Surface 

A11=F11/pU17ra2h 
Solid 

A 
£  =    31 

h  Aii 
Free  Surface 

Boundary  Cond. 

Boundary  Cond. 

Boundary  Con 

10.0 

0.935 

1.006 

0.476 

9.0 

0.908 

1.002 

0.473 

8.0 

0.889 

1.000 

0.469 

7.0 

0.869 

1.000 

0.465 

6.0 

0.844 

1.003 

0.459 

5.0 

0.813 

1.007 

0.453 

4.0  • 

0.769 

1.014 

0.446 

3.0 

0.704 

1.022 

0.437 

2.0 

0.597 

1.031 

0.427 

1.5 

0.515 

1.036 

0.421 

1.0 

0.404 

1.045 

0.416 

0.5 

0.255 

1.069 

0.410 
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TABLE  3 

CYLINDER:  COMPUTER  GENERATED  DATA  (Horizontal  Motion) 
252  Effective  Points 


h 

• 

A11=P11/pU17ra' 
Free  Surface 

"h 

A11=F11/pU1Tra2h 
Solid 

h  AX1 
Free  Surface 

Boundary  Cond 

Boundary  Cond. 

Boundary  Cond 

.0.0 

0.907 

0.999 

0.471 

9.0 

0.893 

0.998 

0.468 

8.0 

0.877 

0.998 

0.464 

7.0 

0.858 

0.999 

0.460 

6.0 

0.835 

1.002 

0.455 

5.0 

0.804 

1.006 

0.449 

4.0 

0.761 

1.011 

0.442 

3.0 

0.697 

1.017 

0.433 

2.0 

0.590 

1.022 

0.424 

1.5 

0.510 

1.026 

0.419 

1.0 

0.400 

1.032 

0.414 

0.5 

0.250 

1.051 

0.410 
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TABLE  6 

CYLINDER:   COMPUTER  GENERATED  DATA  (Rotational  Motion 
about  vertical  axis)  192  Effective  Points 


h 

A66=M66/(pfi37ra2h3) 

A66=M66/(pVa2h3 

Free  Surface 

Solid 

Boundary  Cond. 

Boundary  Cond. 

10.0 

0.266 

0.322 

9.0 

0.258 

0.320 

8.6 

0.248 

0.319 

7.0 

0.237 

0.317 

6.0 

0.225 

0.314 

5.0 

0.209 

0.312 

4.0 

0.190 

0.308 

3.0 

0.166 

0.302 

2.0 

0.132  . 

0.294 

1.5 

0.110 

0.288 

1.0 

0.083 

0.282 
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pl 


TABLE  7 


CYLINDER:   COMPUTER  GENERATED  DATA  (Pressure 

Coefficients,  Free  Surface  Boundary 
Condition) 

Pl 


plula  d  =  d/a  h  =  h/a 


0.929  4.722  5.0 

0.925  4.167  5.0 

0.916  3.611  5.0 

0.900  3.056  5.0 

0.875  2.500  5.0 

0.835  1.944  5.0 

0.769  1.389  5.0 

0.649  0.833  5.0 

0.394  0.278  5.0 


2.833  0.847  3.0 

2.500  0.841  3.0 

2.167  0.827  3.0 

1.833  0.804  3.0 

1.500  0.769  3.0 

1.167  0.716  3.0 

0.833  0.637  3.0 

0.500  0.510  3.0 

0.167  0.279  3.0 
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TABLE   7    cont 


1      P-,U,a  d=   d/a  h=  h/a 


0.743  1.889  2.0 

0.735  1.667  2.0 

0.719  1.444  2.0 

0.693  1.222  2.0 

0.655  1.000  2.0 

0.601  0.778  2.0 

0.522  0.556  2.0 

0.405  0.333  2.0 

0.208  0.111  2.0 


0.522  0.944  1.0 

0.515  0.833  1.0 

0.500  0.722  1.0 

0.477  0.611  1.0 

0.445  0.500  1.0 

0.401  0.389  1.0 

0.340  0.278  1.0 

0.255  0.167  1.0 

0.120  0.056  1.0 
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•      TABLE    11 
SPHERE:      COMPUTER   GENERATED   DATA    (h   =    4) 

d      A^F^pU^I3     A11=F11/pU1jrTa3     A22=F22/ptJ2^ra3     A22=P22/pU2yra3 

Free  Surface   Solid         Free  Surface   Solid 
Boundary  Cond. Boundary  Cond. Boundary  Cond.  Boundary  Cond. 


3.0 

0.663 

0.682 

0.898 

0.904 

2.9 

0.595 

0.613 

0.717 

0.725 

2.7 

0.556 

0.576 

0.620 

0.632 

2.5 

0.538 

0.559 

0.578 

0.595 

2.3 

0.526 

0.551 

0.554 

0.578 

2.0 

0.514 

0.547 

0.529 

0.570 

1.8 

0.507 

0.549 

0.513 

0.573 

1.6 

0.498 

0.555 

0.495 

0.585 

1.4 

0.485 

0.566 

0.470 

0.610 

1.2 

0.465 

0.590 

0.433 

0.665 

1.1 

0.451 

0.613 

0.407 

0.725 

1.0 

0.433 

0.682 

0.375 

0.904 

0.5 

0.190 

0.225 

0.300 

0.880 
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TABLE  12 

SPHERE:   COMPUTER  GENERATED  DATA  (Horizontal  Motion, 
Bottom  Mounted) 


6.0 

5.0 

4.0 

3.0 

2.5 

2.0 

1.90 

1.62 

1.43 

1.22 

1.00 


A11=F11/(PU147Ta3/3)  A11=F11/(pU14Tra3/3) 


Free  Surface 

Solid 

Boundary 

Cond. 

Boundary  Cond. 

0.669 

0.674 

0.667 

0.676 

0.663 

0.682 

0.648 

0.701 

0.625 

0.732 

0.555 

0.901 

0.529 

0.960 

0.440 

0.860 

0.371 

0.752 

0.292 

0.609 

0.211 

0.451 
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TABLE  15 

HEMISPHERE:   COMPUTER  GENERATED  DATA  (Horizontal  Motion) 


A11=P11/(pU12T;a3/3  A11=F11/(pU127Ta3/3) 


6,0 
5.0 
4.0 
3.0 
2.5 
1.7 
1.4 

1.2 
1.0 


Free  Surf 

'ace 

Solid 

Boundary 

Cond. 

Boundary  Cond. 

0.516 

0.518 

0.515 

0.519- 

0.514 

0.520 

0.510 

0.526 

0.505 

0.533 

0.481 

0.567 

0.455 

0.610 

0.421 

0.674 

0.366 

0.901 
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•  TABLE  16 
CONE:   COMPUTER  GENERATED  DATA  (Horizontal  Motion) 


6.0 
4.0 
3.0 
2.0 
1.5 
1.25 


All=pll/pul* 

A11=F11/pU^ 

Free  Surface 

Solid 

Boundary  Cond. 

Boundary  Cond 

0.885 

0.901 

0.869 

0.922 

0.845 

0.960 

0.766 

1.109 

0.664 

1.405 

0.578 

1.857 

6.0 
4.0 
3.0 
2.0 
1.5 
1.25 


Free  Surface 

Solid 

Boundary 

Cond. 

Boundary  Cond 

0.336 

0.336 

0.338 

0.335 

0.339 

0.333 

0.348 

0.324 

0.370 

0.305 

0.394 

0.281 
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APPENDIX  B 
EXPERIMENTAL  DATA 

TABLE  4 
CYLINDER:   CALIBRATION  DATA  (Horizontal  Motion) 

Cylinder  (A) 
Frequency  (HZ)  Added  Weight  (lbs) 

18.00 
16.15 
15.70 
15.30 
14.45 
14.15 
13.85 
13.80 
13.25 
12.70 
11.80 
11.10 


18.95 
16.63 
15.27 
14.86 
14.10 
13.80 
13.48 
12.47 
12.04 

11.67 


0.00 

0.793 

1.022 

1.275 

1.814 

2.043 

2.296 

2.308 

2.790 

3.329 

4.360 

5.388 

Cylinder  (B) 

0.00 

0.864 

1.561 

1.866 

2.351 

2.591 

2.880 

3.888 

4.373 

4.890 
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TABLE  5 
CYLINDER:   EXPERIMENTAL  DATA  (Horizontal  Motion) 


Cylinder  (A) 


Frequency  (HZ)      A,n  = 


Added  Wt . 
11  Displaced  Wt . 


5.30  12.10  0.835 

4. 80  12.50  0.821 

4.20  13.00  0.804 

3.70  13.48  0.775 

3.15  14.10  0.722 

2.75  14.55  0.682 

2.25  15.30  0.618 

2.00  15.66  0.585 

1.65  16.10  0.548 

1.25  16.75  0.441 

0.90  17.20  0.380 

0.40           ■  17.78  0.221 

0.00  18.00  0.000 

Cylinder  (B) 

5.30  12.43  0.818 

4.80  12.80  0.816 

4.20  13.40  0.788 

3".  70  13.92  0.745 

3.15  14.60  0.694 

2.75  15.07  0.678 


68 


TABLE  9 
SPHERE:   CALIBRATION  DATA  (Vertical  Motion) 

Frequency  (HZ)                           Added  Weight  (lbs.) 

19.8  0.0 

19.0  0.3 

18.6  0.5 

18.3  °-7 

17.5  1-0 

16.7  1-3 

16.4  1.5 

16.0  1-7 
15.4  2.0 

15.1  2-3 
15.0  2.5 

14.6  2.7 
14.4  3.0 
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TABLE  10 

SPHERE:  EXPERIMENTAL  RESULTS  (Vertical  Motion,  Infinite  fluid 

Depth) 

.  .     Added  Weight 

22   Displaced  Weight 


6.0  0.512* 

5.0  0.480 

4.0  0.514 

3.0  0.504 

2.0  0.500 

1.75  0.475 

1.5  0.463 

1.25  0.416 

1.0  0.367 
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TABLE  13 
SPHERE:  CALIBRATION  DATA  (Horizontal  Motion,  Bottom  Mounted) 

Frequency  (HZ)  Added  Weight  (lbs.) 

14.65  0.379 

14.35  0.592 

•    14.00  0.828 

13.95  0.864 

13.65  1.077 

13.40  1.313 

13.00  1.625 

12.75  1.873 

12.25  2.390 

11.85  2.875 

11.45  3.397 

10.75  4.405 
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TABLE  14 
SPHERE:  EXPERIMENTAL  DATA  (Horizontal  Motion,  Bottom  Mounted) 

added  wt . 


'11   displaced  wt . 


3.72  12.05  0.644 

3.40  12.00  0.661 

3.13  12.05  0.644 

2.90  12.10  0.624 

2.67  12.15  0.612 

2.50  12.20  0.605 

2.33  12.25  0.585 

2.17  12.35  0.563 

2.00  12.45  0.536 

1.83  12.70  0.465 

I.67  12.90  0.421 

1.50  13.20  0.355 

1.33  13.55  0.287 

1.17  13.85  0.230 

1.00  14.15  0.179 

0.00  15.32  0.00 
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TABLE  17 

CONE:  CALIBRATION  DATA  (Horizontal  Motion) 

Frequency  (Hz)  Added  Weight  (lbs.) 

17.61  0.418 
16.89  0.631 
16.50  0.867 
15.53  1.352 
15.11  1.664 

14.62  1.912 
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TABLE  18 
CONE:  EXPERIMENTAL  DATA  (Horizontal  Motion) 

Frequency  (Hz) 


Added  Wt . 


11   Displaced 
Wt. 


0.79  15.32  0.281 

0.64  16.03  0.214 

0.48  16.89  0.142 

0.33  17.79  0.095 

0.17  18.66  0.065 
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Fig.    13   -   Calibration   curve    ■ 

(Sphere    -  Vertical   Motion) 
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